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Abstract. - 

The two-dimensional Hubbard model with a bimodal distribution of on-site interactions, 
P(C/j) = (1 — /)5(C/j — (7) + /(5(C/j), is studied using finite temperature Quantum Monte Carlo 
and dynamical mean-field theory. Long range antiferromagnetic order off half-filling is stabilized 
by the disorder, due to localization of the dopants on the U = Q sites. Whereas in the clean 
model there is a single gap at n = 1, for nonzero / we find the compressibility and density of 
states exhibit gaps at two separate fillings. 



The Hubbard model exhibits both Mott metal-insulator and magnetic phase transitions. On 
the one hand, the interaction U induces a gap at half-filling by separating many-body states 
with doubly occupied orbitals (upper Hubbard band) from those with holes (lower Hubbard 
band) when U becomes larger than the non-interacting bandwidth W . On the other hand, 
also near half-filhng, there is a tendency to antiferromagnetic (AF) ordering of the electron 
spins. For interacting fermions, the issue has recently been raised whether the Mott transition 
ever occurs in the absence of associated symmetry breaking such as magnetic order Q. 

The introduction of "impurity" sites where U — Q 'm principle could separate AF order 
and the Mott transition. The Mott gap will be shifted to densities greater than n = 1 since 
some sites can be doubly occupied without any on-site repulsion energy cost. Meanwhile, it 
is likely that the Fermi-surface instability responsible for opening the AF (spin density wave) 
gap remains at half-filling. In this paper we will study a model Hamiltonian incorporating this 
effect. Our main conclusions are: {1) U — Q sites can induce long range AF order at densities 
which are spin disordered in the clean model, through localization of the doped particles. (2) 
The dependence of the occupation on the chemical potential and the behavior of the density 
of states exhibit a Mott gap, shifted from half-filling, and also a gap at half-filling resulting 
from induced AF order on the [/ = sites. 

There are a number of experimental systems where the effect of the introduction of nonmag- 
netic impurities has been studied. Examples include doping Zn or Ga for La in La2Cu04, where 
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the critical concentration for the destruction of AF order, Xc ~ 0.10—0.15 is considerably larger 
than for doping with Sr, which is not isovalcnt and also Zn doping in ladder compounds 
where an AF phase is stabilized jsj. Numerical work on the effect of impurities on Neel order 
has focused on spin systems Enhanced local correlations arise from restrictions on the 
singlet bond patterns due to the defects 

We study the electronic lattice Hamiltonian, 

<ij>,(T icr i 

The kinetic energy is described by hopping between nearest-neighbor lattice sites, (ij), on a 
square {N=L x L) lattice. The interaction values C/j are chosen from a bimodal distribution, 
[/j e {0,C/}, with probabilities / and 1 — /, respectively. Upon electron doping the C/ = 
sites will be filled thereby quenching their local moment and eventually acting as magnetically 
inert impurities. 

We use two complementary numerical approaches, finite dimensional "determinant" quan- 
tum Monte Carlo (QMC) algorithm 1^, |^] and the closely related "dynamical mean field theory" 
(DMFT) I, |, 0, also employing QMC 0. The former technique allows an approximation 
free solution of the equilibrium thermodynamic properties, but only on finite lattices and in a 
somewhat restricted range of parameters. The latter technique introduces a local approxima- 
tion to the self-energy, which becomes exact in the limit of infinite dimensions, and allows a 
more complete exploration of parameter space in the thermodynamic limit. DMFT has proven 
to give reliable results for thermodynamic and dynamic properties in d = 2 and d = 3 ^ , 
and has also been applied to study disorder effects on interacting systems [^2[ |l3|. While for 
d = 2 expectation values are averaged over typically tens of (static) disorder configurations, in 
DMFT this average reduces to a weighted sum over the two possible values of Ui, equivalent 
to the coherent -potential approximation fl^ . 

Like the original Hubbard model, Hamiltonian (j^) is particle-hole symmetric at half-filling 
(n = 1 and = 0), i.e. it is invariant under the "staggered" particle-hole transformation 
cI — > (— l)^c. , so there is no minus-sign problem at n = 1 which would preclude simulations 
for large lattices at low temperatures T. Physically, the particle-hole symmetry corresponds 
to different chemical potentials on the two constituents such that at ^ = the local density 
expectation values nj are homogeneous, i.e. independent of J7j |jl^. The choice here of the 
symmetric limit of the model is identical in spirit to the focus in the heavy fermion literature 
on the 'symmetric' Anderson model, in which a similar specific relation between f-site energy 
and f interaction strength is assumed to bring out most clearly the physics of that Hamiltonian. 
In the following we will discuss the effects of disorder on (i) AF order and (ii) the charge gap, 
presenting for each case first d = 2 data which arc then complemented by DMFT results. 



AF order. — First we consider the case n — 1 and study the effect of an increasing 
concentration f oiU = sites on the stability of AF long range order. The static AF structure 
factor S'(7r,7r) is calculated for different lattice sizes (L < 10) and temperatures, and averaged 
over up to 40 disorder configurations. For L < 10, 5(7r, tt) is found to saturate at T « t/8. 
From the saturated values the ground state sublattice magnetization AI can be extrapolated 
using a finite-size scaling according to spin wave theory jl6|, S{n,Tr)/L'^ = M^/3 -I- 0{1/L). 
Scaling plots for different fatU = 8t are shown in the inset of Fig. 1. For / < 0.36, ^(Tr, tt) /L^ 
extrapolates to a finite value in the thermodynamic limit. For / = 0.5 the linear scaling with 
1/L does not hold anymore, since it leads to negative values of AP, indicating the absence 
of AF order in the thermodynamic limit. Fig. 1 presents the extrapolated values of M as a 
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Fig. 1. - Ground state staggered magnetization M as a function of fraction of [/ = sites, /, 
for U = 8t, n = 1. Inset: Finite-size scaling of the antiferromagnetic structure factor. / =0.0, 
0.125, (filled triangles and squares) are ordered. / = 0.500 (open circles) is not. For values of 
/ that do not correspond to an integer number of defects, we have interpolated the results for 
the two bracketing concentrations. 

function of /. A spin wave suppression of the staggered magnetization [l9|| . might be the origin 
of the weak enhancement of M near / = 0.1. 

The U — impurity sites can induce AF order at a doping for which the clean model is 
spin disordered. At finite doping the minus-sign problem prevents the simulations necessary 
for the finite-size analysis. We therefore employ DMFT to calculate the phase diagram in 
the density-disorder plane, shown in Fig. 2 for T = t/8. The phase boundary is obtained 
both from the vanishing of M in the AF phase and from the divergence of the staggered 
susceptibility in the spin disordered (D) phase. The region of AF long range order extends to 
a larger doping for finite /; the presence of impurity sites stabilizes order. In the clean model, 
doping away from half-filling introduces extra particles which are mobile and hence especially 
effective at disturbing correlations over the entire lattice. U — defects provide localizing 
sites which are energetically favorable for these extra particles. Thus, while doping destroys 
moments locally, it has a much smaller effect on long range correlations. In two dimensions, 
the critical doping value for the clean f = model is smaller than in DMFT and probably zero 
because of competing charge or incommensurate magnetic instabilities. However, we expect 
that the enhancement of the AF phase by nonzero / will also be present in d = 2, since the 
underlying mechanism, the defect-induced localization, is observed in d = 2, too. 

Charge gap. — In the extreme strong coupling limit, t = 0, a plot of the density n as a 
function of chemical potential /i exhibits plateaus at n = 1 ± /: n is zero until /i — — C//2, at 
which point n will jump to n = 1 — / (singly occupying sites which have nonzero U.) When 
fi goes through zero, n will jump ton==l — / + 2/ = l + / (singly occupied nonzero U sites 
plus doubly occupied U — sites). Finally, at /i = +U/2, the density jumps to n = 2 (all 
sites doubly occupied). When t is turned on these sharp steps will be rounded by quantum 
fluctuations, but we might expect the compressibility n — dn/dfi to remain large at n = 1 
and remain small at n = 1 + / . In order to explore precisely the effect of nonzero hopping the 
dependence of n on the chemical potential is calculated using d = 2 QMC. Fig. 3 shows results 
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Fig. 2. - Phase diagram in DMFT as a function of filling, n, and fraction of U ~ sites, /. 
For < / < 0.4 the antiferromagnetic (AF) region around half-filling is enhanced relative 
to / = 0. At fixed density n > 1, the introduction oi U — impurities localizes the mobile 
dopants and restores long range order. 



Fig. 3. - Density n vs. chemical potential /i for different values of / on a 6 x 6 lattice with 
U — 8t. Error bars are within the symbol size when not shown [f8]. The Mott plateau 
is shifted away from half-filling. Dashed lines indicate the values 1 + / (corresponding to 
complete occupation oi U = sites). 



for a 6 X 6— lattice at T — t/8 and U — 8t At / = there is a wide charge gap, visible in 
the broad plateau at n = 1 starting at /i = 0. For / > the "Mott" gap occurs close to the 
expected density of n = 1 + / (dashed lines). It terminates at a chemical potential close to 
the fic for / = 0. That is, the chemical potential to force double occupation of the non-zero U 
sites is unaffected by the presence of the U = sites. In addition, however, a second plateau 
remains at half-filling, an effect not predicted by the t = analysis. 

In order to elucidate the nature of this remnant gap at half-filling, we examine the effect 
of the chemical potential on the staggered magnetization M. Fig. 4 shows DMFT results for 
the densities and sublattice magnetizations vs. /j, at / = 0.11 for [/ = and U = 8t sites 
separately The densities show quantitatively the same dependence on /.t in d = 2 and in 
DMFT ||l^, however the sublattice magnetization off half-filling can be obtained in DMFT 
only. The density on the U — 8t sites hardly changes for fi < 1.6t, but then begins to rise, 
at which point M also abruptly vanishes. For small /x there is induced AF order also on the 
U = sites. This AF order leads to a gap even on the U = sites due to the doubling of the 
unit cell. The two different gap values on the U = 8t and U = sites explain the double gap 
structure in the total density (Fig. 3). 

Surprisingly, M becomes negative on the U — sites when n{U — 0) starts to saturate. The 
reason is that an electron on a J7 = site with spin parallel to its neighbors is more strongly 
localized (due to Pauli's principle) than an electron with opposite spin. Thus the net moment 
on the U = site is parallel to its neighbors, i.e. opposite to the total staggered magnetization. 
Within the DMFT approach AF order is more stable against disorder (/^ — 0.75 ± 0.02 at 
T = and n = 1, compared to fc ~ 0.4 in d = 2), as is typical for mean-field theories. 

The charge gap can also be obtained from the density of states (DOS). We employ the 
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Fig. 4. - Density n (solid symbols) and sublattice magnetization M (open symbols) vs. fi in 
DMFT separately for [/ = sites (triangles) and U = 8t sites (squares). For small fi there is 
AF order on both types of sites. For > 1.6t the density on the U — 8t sites increases and 
AF order breaks down. 

Fig. 5. - Density of states in DMFT on [/ = and C/ = 8i sites for T = t/8, / = 0.11 at a) 
half-filling and h) n = 1.11. All spectra have a vanishing DOS at the Fermi energy. 



DMFT approach, combined with the "Maximum Entropy" method , which gives excellent 
agreement with d = 2 and d — 3 results in the clean case ||2^, At half-filling and low T 
(Fig. 5a) the DOS for the U — 8t sites shows the typical AF structure with two broad 
Hubbard bands, two quasi-particle peaks at low energies, and a gap at the Fermi energy, 
Ep = 0. On the U = sites the Hubbard bands appear at lower energies and the gap is 
very small but still present, and will increase at lower T. Thus there is indeed a gap in the 
antiferromagnetic phase, even for the U = sites. For small doping (1.0 < n < 1.11) the 
compressibility was seen to be large (see Figs. 3 and 4). At n = 1.11, where the density on the 
U = sites saturates, the DOS again vanishes at the Fermi energy, Ep — n — lA5t (Fig. 5b). 
At this point there is still AF order. 

Summary. — We can summarize the physics of the charge gap and the AF order in terms 
of the phase diagram of Fig. 2. In the pure case AF order breaks down when the density 
reaches a critical value 7ic(0). This value depends on dimensionality and may very well be 1.0 
in d = 2. At a finite concentration /, the defect sites localize dopants, thereby keeping the 
rest of lattice at half filling, until they are (almost) doubly occupied. For small concentration 
/ the defect sites hence stabilize AF order up to a total density of ndf) ~ nc(0)(l — /) + 2/. 
Eventually, at large enough /, rif, is driven to one and AF order ceases to exist even for n — I. 
The density exhibits a plateau associated with AF order, which exists on both the nonzero 
U and the U = sites, as the chemical potential is changed across half-filling (n = 1); and 
a second "Mott" plateau, associated with double occupancy, as one crosses the high filling 



6 



EUROPHYSICS LETTERS 



boundary of the AF region. 
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